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Abstract 

The production of the Xc2 charmonium state in two-photon collisions has been 
measured with the Belle detector at the KEKB e + e~ collider. A clear signal for 
Xc2 —> 7 J J/tp — ► is observed in a 32.6 fb _1 data sample accumulated at 

center-of-mass energies near 10.6 GeV, and the product of its two-photon decay 
width and branching fraction is determined to be r 77 (x C 2)<6(Xc2 - ► 'yJ/ip)B(J/'ip — > 
£+£-) = 13.5 ± 1.3(s*at.) ± l.l(syst.) eV. 

Key words: two-photon collisions, charmonium, Xc2, partial decay width 
PACS: 13.20.Gd 13.40.Hq 13.65,+i 14.40.Gx 
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1 Introduction 



The two-photon decay widths (T 77 ) of the even charge-parity charmonium 
states provide valuable information for testing models describing the nature 
of heavy quarkonia. Various theoretical calculations describing the quark- 
antiquark system predict the value of r 77 (x C 2) to be within the range 0.2 
- 0.8 keV [1]. A precise experimental determination of r 77 (x c2 ) will provide 
a strong constraint on these models. The ratio of the two-gluon decay width 
to the two-photon decay width r 99 (x c2 )/r 77 (x c2 ) has been calculated within 
the framework of perturbative QCD with first-order correction [2] and the 
result gives r 77 (x c2 ) = 0.47 ± 0.13 keV.Q A comparison of r 77 (x c2 ) with the 
two-gluon width will provide a way to study the validity of perturbative QCD 
corrections for quarkonium decays. 

Measurements of two-photon decay widths for charmonium states are difficult 
because of their small production cross section and small detection efficien- 
cies. To date, several experiments have reported the observation of two-photon 
production of the Xc2 in the decay channel Xc2 —> lJ/4>, J 1 4 s — ► where 
I = e or /i [5]. This channel is suitable for the experimental determination of 
1^77 {Xco)i since the decay branching fraction is known with a relatively small 
error, B(xc2 -> jJ/i/j)B(J/ip -> t+£~) = (1.59 ± 0.13)% for both lepton fami- 
lies [3]. However, the two-photon decay width results obtained from previous 
measurements with this process seem to be systematically larger than those 
frompp — > Xc2 77 experiments [6]. Further studies with high statistics data 
samples are needed to clarify the situation. 

Recently, the CLEO collaboration has reported a measurement of r 77 (x c2 ) in 
the 7T + 7r~7r + 7r~ final state [7]. Although the branching fraction of this decay 
mode, B(xc2 — ► 7r + 7r~7r + 7r~) = (1.2 ± 0.5)% [3], is comparable to that of the 
jJ/ip — > £ + £~ , y mode, its large error precludes a precise determination of 

r 77 (xc2)- 

We have measured Xc2 production in two-photon processes using the decay 
channel Xc2 — * 7^/V? J/*l> ~ > ■ The results are based on a 32.6 fb _1 data 
sample collected with the Belle detector. 



1 This value is derived from measurements [3] using the assumption that T gg (xc2) 
is given by T(x c2 — > hadrons) — T(xd ~^ hadrons) [4]. Here the strong coupling 
constant is set to be a s (m c ) = 0.29 ± 0.02. 
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2 Experimental data and detector system 

The experiment was performed with the Belle detector [8] at the asymmetric 
e + e~ collider KEKB, where an 8.0 GeV e~ beam collides with a 3.5 GeV e + 
beam with a crossing angle of 22 mrad. We use a 29.6 fb _1 sample of data 
collected at the cm. energy corresponding to the peak of the T(4S) resonance 
(10.58 GeV) and a 3.0 fb _1 sample collected 60 MeV below the peak. 

The basic topology of events that we select is two tracks of opposite charge and 
a photon. The recoiling e + and e~ are not tagged in order to select quasi-real 
two-photon collisions with high efficiency Events induced by highly virtual 
photons (i.e., photons with high-Q 2 ) are effectively rejected by a strict trans- 
verse moment um(p t ) requirement applied to the Xc2 daughter particles, as 
described in the following section^ 

The charged track momenta are measured with a cylindrical drift chamber 
(CDC) located in a uniform 1.5 T magnetic field. Track trajectory coordinates 
near the collision point are provided by a silicon vertex detector (SVD). Photon 
detection and energy measurements are performed with a Csl electromagnetic 
calorimeter (ECL). The resolutions of track momentum and photon energy 
measurements are 0.4% for the leptons, which, for the signal process, have 
typical p t values of 1.4 GeV/c, and 2.0% for the photons, with typical energies 
of 0.4 GeV. The magnet return iron is instrumented to form the Kl and muon 
detector (KLM), which detects muon tracks and provides trigger signals. 

The majority of the signal events are triggered by two-track triggers that 
require at least two tracks with transverse momenta larger than 0.2 GeV/c 
detected in the CDC in coincidence with matching signals from TOF counters 
and trigger scintillation counters, isolated cluster or energy-sum signals from 
the ECL, or muon tracks in the KLM. A constraint on the opening angle in the 
plane transverse to the e + beam axis {rip plane), Aip > 135°, is applied at the 
trigger level. The signal inefficiency due to this trigger constraint is negligibly 
small because the final-state leptons tend to be back-to-back in the rip plane 
due to the kinematic properties of our selected final states, namely the strict 
Pt requirement on the Xc2 and the small mass difference between the Xc2 and 
J/ip. The events from the J/ijj — > e + e~ mode are efficiently triggered by a 
total energy trigger derived from the ECL with a threshold set at 1.0 GeV. 

The main backgrounds are leptonic final states from QED processes such as 
e + e~ — > e + e~£ + £~ , y. The E/p information, which is the ratio of the energy 
deposit in the ECL to the track's momentum, is used to identify the leptons 

2 Q 2 is defined as the negative of the invariant mass squared of a virtual incident 
photon. It is approximately equal to \pt\ 2 of the virtual photon with respect to the 
e + e~ beam axis. 
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and eliminate small backgrounds that contain hadron tracks. 



3 Event selection 

The event selection criteria are as follows: (1) Exactly two oppositely charged 
tracks reconstructed by the CDC, where both tracks satisfy the following 
laboratory frame conditions: —0.47 < cos# < +0.82, where 9 is the polar 
angle; p t > 0.4 GeV/c; \dr\ < 1 cm, \dz\ < 3 cm, where (\dr\, dz) are the 
cylindrical coordinates of the track's point of closest approach to the nominal 
collision point in the rip plane; \Adz\ < 1 cm, where Adz is the difference 
between the dz's of the two tracks; and no other well reconstructed tracks 
with p t higher than 0.1 GeV/c. (2) The opening angle (a) of the two tracks 
satisfies cos a > —0.997. (3) There is just one electromagnetic cluster in the 
ECL with an energy E 1 > 0.2 GeV and isolated from the nearest charged track 
by an angle greater than 18°. (4) The scalar sum of the momenta of the two 
charged tracks is less than 6 GeV/c, and the invariant mass of the two tracks 
is between 1.5 and 4.5 GeV/c 2 . (5) The total energy deposited in the ECL is 
less than 6 GeV. (6) The absolute value of the total transverse momentum 
vector in the cm. frame of the e + e~ beams, |p* tot | = \pt + + Pt + P* 7 |; 
is less than 0.15 GeV/c, while that for the two tracks only, |p^ + + pl~\, is 
larger than 0.10 GeV/c, where pl + , pl~ and \\p are measured transverse 
momentum vectors (defined as two-dimensional momentum vectors projected 
onto the plane perpendicular to the beam axis in the e + e~ cm. system) for 
the positive track, the negative track and the photon, respectively. (7) For 
electron pairs, both tracks are required to have E/p > 0.8; for muon pairs, 
both tracks are required to have E/p < 0.4. 

Selection criterion (2) rejects cosmic-ray backgrounds. Criterion (6) rejects 
two-photonic lepton-pair production events with radiation from a recoil elec- 
tron or with a fake photon, which tend to populate the region |p^ + + p^~| ~ 0. 
The scatter plot in the |p^ + + F*t _ I - 1 Pt tot I plane before the application of re- 
quirement (6) is shown in Fig. 1(a), where two separate clusters of events at 
|Pt + + Pt~| ~ and |P( tot | ~ are apparent. The latter cluster corresponds to 
exclusive final states produced by two-photon collisions. The distribu- 

tion from the signal events obtained from the Monte Carlo (MC) simulation 
(described in Sect. 5) is shown in Fig. 1(b). 

We correct the absolute momenta of detected electrons or positrons for 
bremsstrahlung in e + e~7 event candidates. If photons of energy between 0.02 
and 0.2 GeV are present within a cone of half-angle 3° around the electron 
direction, the energy of the most energetic photon in the cone is added to the 
absolute momentum of the track. This correction is also effective for partially 
compensating for the radiative-decay events of J/ip, J/ip — > e + e~^. 
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(b) 



0.2 0.4 0.6 0.8 1 

|p{++pf-|(GeV/c) 

Fig. 1. Scatter plots for the absolute values of the two kinds of vector sums of the 
transverse momenta: the horizontal axis for the sum of the two charged particles, 
and the vertical axis of all three particles including the photon: (a) for real data and 
(b) for Monte Carlo events of the signal process. Straight lines show the selection 
requirements. 

4 Derivation of the number of signal events 

Figure 2 shows a scatter plot of the invariant mass of the two tracks (M H ) ver- 
sus the invariant-mass difference AM = M + _ 7 — M + _ for the selected events, 
where M + „ 7 is the invariant mass of all three particles. A clear concentration 
is observed around M + _ = 3.097 GeV/c 2 and AM = 0.459 GeV/c 2 , the signal 
region for \c2 — > iJ/^- 

The mass difference distribution is shown in Fig. 3(a) for the events falling 
within the J/ip signal mass region 3.06 < M + _ < 3.13 GeV/c 2 . After the final 
selection requirement for signal candidate events, 0.42 < AM < 0.49 GeV/c 2 , 
is applied, 176 events remain. Of these, 82 events have electron pairs and 94 
have muon pairs. The contribution of Xci production, which would peak at 
AM = 0.42 GeV/c 2 , is estimated to be less than one event, as expected from 
the suppression of spin-1 meson production in quasi-real two-photon collisions 
[9]. 
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Fig. 2. Scatter plot of the invariant mass M_| of the two-track system versus the 

mass difference AM = 7 — The arrows indicate the nominal J/ip mass 

and the AM value expected for \c2 - ► lJ/^- 

The AM distribution of events in the J/^-mass sideband regions (2.65 < 
M+_ < 3.00 GeV/c 2 and 3.15 < M + _ < 3.50 GeV/c 2 ) is used to determine the 
background contribution (Fig. 3(b)). The AM distribution for J/ip sideband 
events between AM = 0.15 and 0.70 GeV/c 2 is fitted by an appropriate 
function^] to obtain the shape of the background distribution. We normalize 
this to the AM distribution of events in the J / ip signal region, and determine 
the number of background events in the Xc2 signal region to be 40.0 ± 2.7 
events. This normalization fit uses only data in the Xc2 sideband regions, 
0.15 < AM < 0.29 GeV/c 2 and 0.50 < AM < 0.70 GeV/c 2 . Since we expect 
9±7 events around AM = 0.32 GeV/c 2 from the Xco [3,7], we avoid that mass 
region. This expected Xco yield is consistent with the small excess of events 
seen near AM = 0.32 GeV/c 2 in Fig. 3(a). We determine the number of Xc2 
signal events to be 136.0 ± 13.3 after subtracting the number of background 
events from the total in the signal region. 

We find that the events in AM sideband regions are dominated by non-J/ip 
backgrounds since the quantity and shape of the distribution agree with that of 
the J /if) sideband. The 368 events in the 0.42 < AM < 0.49 GeV/c 2 region in 
Fig. 3(b) consist of 176 electron pairs and 192 muon pairs. These backgrounds 

3 We used the following empirical function for the fit of the background shape: 
A(AM — a)~ b /(l + e -c ( AM_( ')) 1 where a, b, c and d are free parameters, and A the 
normalization parameter. We have confirmed that the change of the AA1 distribu- 
tion shape is small at different M^ points in the sideband region. 
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0.6 0.7 
AM(GeV/c2) 



Fig. 3. The mass difference distributions for (a) events in the J/^-mass region 
(closed circles with error bars) and (b) sideband events. The curves in (a) and (b) 
indicate the results of the fits that are used to determine the background contribu- 
tion in the signal region. The histogram in (a) shows the distribution of the signal 
MC events normalized to the observed signal. The arrows show the signal region 
(a) and the background control region (b). 



are consistent with higher-order QED events such as e + e~ — > e + e~£ + £~7, 
which would give comparable numbers of events with electron and muon pairs. 
In contrast, hadron production (in which pions would fake muons) would give 
primarily events containing muon pairs. Since a complete calculation of this 
process that takes interference effects into account is not available, we cannot 
estimate the background yield theoretically. 

Figure 4 shows distributions for A<p, |p* tot |, | cos^ -7 ! and cos 8^ for the fi- 
nal candidate events, where Atp is the azimuthal angle difference between the 
two lepton's momenta in the laboratory frame, and |P( tot | is the transverse 
momentum of the £ + C.~ r j system in the cm. frame of the e + e~ beams. #+~ 7 
and are the polar angles of the photon in the £ + £~7 cm. frame and of 
the negatively charged lepton in the £ + £~ cm. frame, respectively, where the 
polar angles are measured with respect to the incident e~ direction. The ex- 
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Fig. 4. Comparison of the final samples (closed circles with error bars, backgrounds 
included) with the sum of the signal MC events (open histogram) and estimated 
background contributions (hatched region) for: (a) Aip; (b) |Pf tot |; (c) |cos#+~ 7 j; 
and (d) cos#l~. There are no entries in the Ay? < 135° region in (a). The dashed 
histograms in (c) and (d) show the distribution for the pure helicity=0 production 
case (see the text in Sect. 5). 



perimental distributions are compared with the sum of the signal MC events, 
described in the next section, and the expected background contributions de- 
termined from the J/ ip sideband events, normalized to the observed numbers 
of events. 

It is apparent from Fig. 4(a) that the final-state leptons have a back-to-back 
topology in the np plane, and the experimental data are consistent with the 
MC distribution. The data in Fig. 4(b) show a peak at very small pt values 
( I F** tot I < 30 MeV/c) of the final-state system that is a typical feature of exclu- 
sive production in two-photon collisions; the data are in good agreement with 
the MC prediction. The distribution for sideband events has a less prominent 
concentration near |pj tot | = 0. In the angular distributions of Figs. 4(c) and 
4(d), the data are consistent with the MC. 
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5 Monte Carlo calculations 



We used Monte Carlo (MC) simulated e + e~ — > e + e~Xc2, Xc2 — ► J/4> — ¥ 

£+£- (£ = e or //) events to calculate the efficiency for the signal process. The 
TREPS MC program [10] is used for the event generation. The effects of J/ip 
radiative decays are modeled with the PHOTOS [11] simulation code, which 
generates photon radiation from a final-state lepton generated by TREPS with 
a probability determined by a QED calculation. All of the final-state particles 
in the MC events are processed by the full detector simulation program. 

We assume that the Xc2 decay to the 'yJ/ip final state is an El transition, 
since experimental observations indicate that this transition dominates the 
decay [12]. Since the helicity state of Xc2 produced in two-photon collisions 
is not known, we assume a pure A = 2 state [13], where A is the helicity of 
Xc2 with respect to the 77-incident axis. The measurement of the polar-angle 
distribution of the photon shown in Fig. 4(c) can be used to evaluate the 
possible contribution from a A = component. The A = 2 component produces 
a cos#+ -7 distribution that is proportional to [1 + cos 2 #+ -7 ], whereas the 
A = component is proportional to [5 — 3cos 2 #+~ 7 ]. The dashed histograms 
in Figs. 4(c) and 4(d) show the expected distributions from the pure A = 
state. The present experimental cos#+~ 7 distribution has better agreement 
with a pure A = 2 hypothesis {x 2 /dof = 12.0/9) than that for pure A = 
(x 2 /dof = 45.0/9), or any other mixture of the two helicity states. Here, each 
X 2 /dof is calculated from the |cos#+~ 7 | distribution divided into ten bins, 
where the total event number in the expected signal distribution is fixed to 
the number of observed signal events. 

The trigger efficiency is experimentally determined using Bhabha and fi + fi~ 
events that are collected with two or more different redundant triggers, as 
described above in Sect. 2. We estimate the probability for the signal-process 
events that survive the selection criteria to pass the trigger conditions to be 
(99 ±1)% ((94 ± 3)%) for ee-f (^7) events; and (96 ± 2)% in average. This 
estimation agrees with the results of a trigger simulation that is applied to the 
signal MC events. The trigger efficiency is also confirmed by the experimental 
yield of two-photonic lepton-pair events, e + e _ — > e + e~£ + £~ , which agrees with 
the expectation from a QED calculation [14]. 

The ECL photon energy resolution was studied by comparing experimental 
and MC mass-difference distributions for D*° — > _D°7 decays in e + e~ annihi- 
lation events. We find that the photon energy resolution is 1.3 times the MC 
prediction. The same tendency is confirmed in the Xci — * samples from 

.B-meson decays and rj — > 77 samples from two-photon collisions. If we use 
a correspondingly wider AM distribution for the Xc2 — ► 1-J r /ip, it decreases 
the efficiency of the AM selection by 3.8% from the MC-determined value. 
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We take this effect into account as a correction for the efficiency and assign 
a systematic error of the same size (±3.8%). We confirm that the observed 
position and width of the signal peak in the data are consistent with their 
expected values. 

When the trigger effects and the difference between the data and MC photon 
energy resolutions are taken into account, the overall efficiency is found to 
be 6.6%. This is the average of the e + e~ and yU + /i" decay channels; the ratio 
of the efficiencies for the two lepton species is ee/pL/i = 0.70 ± 0.03. The 
lower efficiency for e + e~ is due to the occasional presence of extra high-energy 
photons from radiative J/ip decays and electron bremsstrahlung. When the 
ee/nn ratio for the background component determined from the J/ip sideband 
events is taken into account, we expect the event yields in the final samples to 
have a ratio ee/pL/j, = 0.74 ± 0.04. This value is consistent with the observed 
ratio for the final experimental samples: 0.87 ± 0.13. 

The two-photon luminosity function is also calculated by the TREPS program 
[10]. For consistency, the same upper cutoff value of the photon Q 2 , Q^ ax = 
1.0 GeV 2 , and the same vector-meson pole effect are used in the calculation 
of the luminosity function and in the event generation. The uncertainty in the 
luminosity function due to the vector- meson pole effect (we adopt the J/tp 
mass) is small, about 2%, since we apply a strict | p>^ tot | requirement that is 
well below the mass scale of the vector mesons. The Q^ax value does not affect 
the product of the luminosity function and the detection efficiency, since it is 
chosen to be large enough to cover the acceptance of our | tot I requirement. 
The value of the TREPS luminosity function is compared with that obtained 
from a full-diagram calculation of the e + e~ — > e + e~fi + fi~ process [14] in the 
small |P( tot | region. From the difference of the two results, the systematic error 
for the luminosity function is estimated to be 5%, which includes ambiguities 
from the choice of the form factor and the finite | p>* tot | requirement. 



6 Results and discussion 

The two-photon decay width of the Xc2 is related to the signal event yield as 



where / Cdt is the integrated luminosity, r] is the efficiency, m Xc2 (=3.556 GeV/c 2 ) 
is the Xc2 mass and L 77 (m Xc2 )(= 7.75 x 10~ 4 GeV" 1 ) is the two-photon lumi- 
nosity function at the Xc2 mass. The total width of Xc2 (2.00 ± 0.18 MeV [3]) 



Yield 
J Cdt 




r 77 (x c2 )£(x c2 - 7^)S(^ - 



1 j/^)B(j/^^£ + n, 
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is much smaller than the present AM resolution (~ 9 MeV), and does not 
affect the present measurement. 

The observed number of events, 136.0 ± 13.3, implies the result r ' 77 (Xc2)£>(xc2 
-> -fJ/tp)B{J/ip -> = 13.5 ± 1.3 ± 1.1 eV, where the first and sec- 

ond errors are statistical and systematic, respectively. This result corresponds 
to r 77 (x c2 Mx c2 -> 7^) = H4 ± ll(stat) ± 9(syst.) ± 2(B.i2.) eV or 
r 77 (x C 2) = 0. 85 ±0. 08(stat. )±0. 07 (sy st. )±0. 07 (B.R.) keV, where the last er- 
rors correspond to the uncertainties of the branching ratios, B(xc2 lJ/i>) = 
(13.5 ± 1.1)% and B{J/^ -> £ + £~) = (11.81 ± 0.20)% [3]. The systematic er- 
ror has contributions from the trigger efficiency (2%), lepton identification 
efficiency (1.5%), photon detection efficiency (2%), inefficiency due to fake 
photons (less than 2%), J/ip detection efficiency (2%), the AM cut efficiency 
(3.8%), background subtraction (2.3%), the luminosity function (5%) and 
other sources (less than 3%); these total 8% when combined in quadrature. 

The error in the background subtraction is derived from the difference in 
signal yields in the AM signal region between the present method (counting 
the events in the signal region and subtracting the background contribution) 
and an alternative method in which the signal and the background components 
are simultaneously fitted to the AM distribution with all the shape and size 
parameters for the background and signal distributions allowed to float, with 
the Crystal-Ball line shape [15] used for the signal distribution. The error of 
the background normalization is also combined with this error. 

The inefficiency due to an extra (fake) photon with E > 0.2 GeV is estimated 
to be less than 2% from an experimental study of the pj-balanced e + e~ — > 
e + e - /i + yU~ process. 

The uncertainty due to the assumption that the Xc2's are produced in a pure 
A = 2 state and decay via pure El transitions is not included in the systematic 
error. Production in the A = state at the 10% level would increase the 
detection efficiency by 7% and decrease the measured r 77 (x c2 ) value by the 
same amount. Meanwhile, a small mixture of the M2 transitions as has been 
indicated by a measurement, M2/E1 ~ — 0.093lo;oii(= «2(Xc2)) in amplitude 
[12], gives only a 2% effect on the efficiency. 

The Belle result for T 77 (x C 2) is compared with those from previous experiments 
[5-7] in Fig. 5. The present result has the smallest statistical and systematic 
errors of all the two-photon measurements and is consistent with the previous 
two-photon results. However, it is larger than the pp results. A review of the 
experimental results of various branching ratios of ip{2S) and Xc decays [16] 
suggests that this discrepancy may come from incorrect values of B(xc2 — ► 
7J/VO and B(xc2 —> Pp) that are used for the derivation of r 77 (x c2 ) in these 
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Fig. 5. Comparison of the Belle result for T 77 (x C 2) value with those from previous 
measurements [5— 7]. The solid error bars show the statistical errors. The length of 
the dashed part in each error bar corresponds to the size of the systematic error, 
including the branching ratio uncertainty. 

experiments. 



7 Conclusion 

We have measured \ci production from two-photon collisions with a 32.6 fb _1 
data sample collected with the Belle detector at the KEKB e + e~ collider, using 
the decay mode Xc2 - > iJ/i 3 , J 1 4 s — > We find 136.0 ± 13.3 signal events 

after background subtraction. The observed polar-angle distributions of the 
photon and leptons are consistent with those expected from the production of 
Xc2 in the pure helicity 2 state. The product of the two-photon decay width of 
Xc2 and branching fractions, T 11 (xc2)B( y Xc2 —> 'J r / , 4>)B(J TV — * = 13.5 ± 

1.3(stat.) ± l.l(syst.) eV, is obtained. This result corresponds to r 77 (x c2 ) = 
0.85 ± 0.08(stat.) ± 0.07(syst.) ± 0.07{B.R.) keV, where the product of the 
branching fractions, B(xc2 -> lJ/4>)B{J/4) -> t+f) = (1.59 ± 0.13)% [3] is 
used. 
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